Unexpected decoupling of stretching and bending modes in protein gels 
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We show that gels formed by arrested spinodal decomposition of protein solutions exhibit elastic 
properties in two distinct frequency domains, both elastic moduli exhibiting a remarkably strong 
dependence on volume fraction. Considering the large difference between the protein size and the 
characteristic length of the network we model the gels as porous media and show that the high 
and low frequency elastic moduli can be respectively attributed to stretching and bending modes. 
The unexpected decoupling of the two modes in the frequency domain is attributed to the length 
scale involved: while stretching mainly relates to the relative displacement of two particles, bending 
involves the deformation of a strand with a thickness of the order of a thousand particle diameters. 

PACS numbers: 83.80.Kn, 82.70.Gg. Published in: Phys. Rev. Lett. 110, 058303 (2013). 



One of the overarching goals in soft condensed mat- 
ter physics is to establish a rigorous connection between 
structural and mechanical properties of a material, and 
to relate those properties to control parameters such as 
the volume fraction, 0, and the interaction strength, U, 
in a colloidal suspension. Particular attention has been 
given to the dynamical arrest transition from a fluid to 
a solid-like state [1]. It has been shown that weakly at- 
tractive colloids {U ~ kBT) at intermediate values of 
(j> (~ 0.1 — 0.4) can form stress-bearing networks as a 
result of the interplay between spinodal decomposition 
and structural arrest pl-fT3]. This represents an uncon- 
ventional route to gelation, where the arrest occurs when 
the dense domains formed during spinodal decomposi- 
tion undergo structural arrest and phase separation can 
no longer proceed, the system being arrested as a self sup- 
porting network. In spite of the interest towards spinodal 
networks for many applications, there is still very little 
understanding of their mechanical behavior. An exam- 
ple of such systems is the solution of the globular pro- 
tein lysozyme [141 , where the effective attraction between 
the proteins depends on temperature, which allows one 
to continuously explore the transition from homogeneous 
liquid to phase separated and arrested states [15]. This 
makes lysozyme a convenient system to investigate the 
interplay between structural evolution and arrest mech- 
anisms. Previous studies showed that the gel states dis- 
play correlation lengths much larger than the particle 
size [IS] but how such feature determines the mechanics 
remained unclear. 

In this letter, we show that the rheological response 
of arrested lysozyme systems displays two distinct elas- 
tic plateaus that are well separated in time. Both mod- 
uli exhibit a dramatic dependence on volume fraction. 
Starting from these observations, we formulate a simple 
model where the two distinct moduli arise from stretch- 
ing and bending modes of the coarse network structure, 
consisting of strands that are significantly bigger than 



the constituent particles. Our results provide guidance 
for the design of amorphous colloidal systems with elastic 
moduli that are remarkably sensitive to changes in 0. 

Our samples consist of hen egg white lysozyme (Fluka, 
L7651) with a radius r ~ 1.7 nm that are suspended in 
20 mM HEPES buffer at pH =7.8, where lysozyme carries 
a net charge of -|-8e. To screen the electrostatic repul- 
sion we add sah, [NaCl]=500 mM [5]. The total volume 
fraction cf) is obtained from the protein concentrations c 
measured by UV absorption spectroscopy using 4> = c/p, 
where p=1.351 g/cm^ is the protein density. As shown 
in previous work the phase behavior of our systems is de- 
termined by temperature [SI I15j : the liquid- liquid phase 
boundary found at lower temperatures is metastable with 
respect to the liquid-crystal coexistence curve observed 
at higher temperatures. The spinodal decomposition is 
found to be arrested for quenches below the arrest tie line 
at r=15 °C, which leads to a bicontinuous state with 
a protein-poor fluid network and a protein-dense glassy 
network, whose features depend on </> and T 

We measure the rheological properties of such glassy 
networks with a commercial rheometer using a cone and 
plate geometry where the temperature is controlled by a 
Peltier-element. In a typical experiment we load a sam- 
ple in the rheometer at a temperature well above the 
binodal line at T = 25 °C. The temperature is then low- 
ered to the final quench temperature at T=10 °C, below 
the arrest tie line. The samples are left to rest for 300 s 
during which the structures stabilize [15j . This procedure 
enables us to reach reproducible steady states with negli- 
gible aging. We explore different arrested states by vary- 
ing the initial volume fraction 0; any experiment where 
crystals form is disregarded. The viscoelastic properties 
of the arrested states are characterized in the linear re- 
sponse regime using two rheology experiments: (i) oscil- 
latory strain experiments yielding the elastic modulus G" 
and the loss modulus G" as a function of the frequency 
/, and (ii) creep tests yielding the compliance J, as a 
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function of time t, in response to a step stress. 

The results of both experiments compare well with 
each other, as shown for (j)= 0.11 after a quench to 
T = 10 °C in Fig 1; the full squares and the open cir- 
cles denote respectively G' and G" obtained in oscillatory 
strain experiment, while the lines correspond to G" and 
G" obtained by converting the data obtained in creep 
experiments [161 . Strikingly, the mechanical response 
function of our system is characterized by two frequency 
domains where the elastic component G' becomes dom- 
inant; these are marked by Goo ^nd Gq at respectively 
the high and low frequency end of the spectrum shown 
in Fig. 1. Both characteristic elasticities are strong func- 
tions of the volume fraction, as shown in Fig. 2, where we 
display representative examples of the spectra obtained 
for different (p in creep experiments. Varying </> by only 
a factor of three leads to a variation of Goo of four or- 
ders of magnitudes. To determine both characteristic 
elasticities with reasonable accuracy we use the scaling 
procedure described in EPAPS [T7] and report the re- 
sults obtained in Fig. 3. Surprisingly, Goo and Go in- 
crease nearly exponentially with (/), where Goo increases 
more quickly with (j) than Gq, such that the two mod- 
uli appear to reach the same magnitude at large enough 
volume fractions. Such exponential dependence is atypi- 
cal for colloidal gels, where the ^-dependence is generally 
well described by power-laws [HI El [T5H^. By con- 
trast, stronger than power law dependencies of the elastic 
moduli on (f> have been observed in ceramics and other 
porous media [21 . 

To rationalize our findings we thus consider the pos- 
sibility to describe our systems as porous media. The 
structure of our systems has been probed in various scat- 
tering experiments [151 122j , from which we can infer that 
the arrested states are coarse network structures that 
are characterized by a correlation length of the order of 
^ ~ 2.5/im. At somewhat smaller length-scales, larger 
scattering vectors q, we find the hallmarks of the Porod- 
regime, consistent with the gel being a bicontinuous net- 
work with a sharp interface between the dilute phase and 
the dense glassy phase. For q > 0.15nm~^, S{q) fol- 
lows a different regime, compatible with the presence of 
heterogeneities at small length-scales. Clearly, while the 
coarse network is characterized by a length scale that 
is ~1000 times larger than the particle size the internal 
heterogeneities extend only over a few particle diameters. 
This large difference in size indicates that the system may 
be considered as almost homogeneous within the strand, 
being mainly characterized by large length scale hetero- 
geneities, which is reminiscent of porous media. 

We then may conceive that the large mismatch be- 
tween length scales is also at the origin of the unusual 
elastic response observed. Let us consider that the me- 
chanical deformation results in an elastic response where 
both stretching and bending modes are involved |23| . In 
principle, the response should be dominated by the softer 
modes, which are the bending modes. However, unbend- 
ing a strand requires the cooperative displacement of sev- 
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FIG. 1. Frequency dependence of the storage G' and loss 
modulus G" obtained for a lysozyme system with (f> = 0.11 
after a quench to T = 10 °C. Symbols denote the data ob- 
tained in oscillatory strain experiments, G' (■) and G" (o). 
Lines correspond to G' and G" obtained by the conversion of 
the inverse compliance measured in a creep experiments 
under identical conditions. (t) is displayed in Fig. |2| 
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FIG. 2. Representive examples of the time evolution of the in- 
verse compliance J{t)^^ obtained in creep experiments using 
stresses that are well within the linear regime. The responses 
are measured after a quench to r=10 °C for (f} = 0.171 (■), 
0.149 (•), 0.13 (A), 0.11 (T), 0.093 (♦) and 0.075 (-<). 



eral hundreds of particles, which may be delayed as com- 
pared to stretching modes which only involve the relative 
displacement of two particles. It is thus conceivable that 
we probe stretching modes at high /, while probing the 
bending modes at low /. 

To test this assumption we consider that the stretch- 
ing modulus Gs is the main contribution to the high 
frequency response Goo. We start from Gm the shear 
modulus of a homogeneous assembly of proteins and es- 
timate the correction due to the presence of pores. Given 
the lengthscale separation between the protein size (nm) 
and the strand size (/im) we can use an effective medium 
theory to account for the effect of the pores on the elas- 
tic properties [24j . Starting from the stress contribution 



3 



of a single inclusion to a homogeneous material, the ef- 
fective medium theory allows us to calculate the elastic 
constants of a porous medium in analogy to Einstein's 
method for calculating the viscosity of a dilute suspen- 
sion of spheres [53] and later extended by S. Arrhenius 
to higher (f) |26j . Because of the assumption of dilute and 
non-interacting pores, one first obtains a linear depen- 
dence of the properties of the porous material on those 
in absence of porosity: Gg — Gmi^ ~ 9^), with e the vol- 
ume fraction occupied by the pores and g a geometric pre- 
factor related in general to the pore geometry; it can vary 
from g = 5/3 for spherical pores up to 15 for ellipsoidal 
pores upon increasing the ellipsoid's anisotropy pill?/] . 
A small increment in the porosity leads to an increment 
of Gs- dGs = —gGgde. Upon integration with the condi- 
tion Gs — Gm at e = 0, one obtains Gg — Gm&~^^ [H]- 
The porosity e is obtained from the volume fraction oc- 
cupied by the dense network phase h, e = 1 — h. In our 
system h increases linearly with <j), as shown in Fig 4a 
where h ~ acf) and a ~ 3.5 [6]. We thus expect that Gg 
scales as 



Gs 



(1) 



For the elastic response at the low-frequency end of 
the spectrum we assume that elasticity is dominated by 
the strand elasticity and hence, because of the strand 
anisotropy, by the bending modes [23j [28] . If we consider 
the strand as an elastic rod of length ^ and diameter m 
[55], its bending energy is K^x"^, where x is the bending 
displacement and the bending constant. The bending 
energy can be also expressed in term of the bending mo- 
ment A, i.e. the moment of forces due to internal stresses 
on a given cross-section of the rod. For small, pure bend- 
ing deformations, A = {E / R){-k / A){m/2Y [29^, where E 
is the Young modulus and R the curvature of the bent 
rod, that is i? ~ ^. Moreover, by schematizing the strand 
as a chain of connected springs with an overall radius of 
gyration ~ ^, the bending displacement due to a force 
applied to the end of the chain is also of the order of 
X ~ ^ [30] ■ As a consequence, the bending modulus is 
Gb ~ Kb/^ ~ A/^^ ~ {m/^)^E. Using the shear modu- 
lus of the strand Gs and eqjljone obtains 



Gb c {m/O^Gs (m/0"e^«^. 



(2) 



For our system ^ is almost independent of 0, as shown 
in Fig 4b. However, since h linearly depends on (j) the 
strand thickness m must depend on cj). Considering that 
the volume fraction occupied by one strand is V ~ ^m^ 
and assuming that the number density of strands is 1/^'^, 
from the volume fraction occupied by the dense phase 

~ V/^^ we obtain the thickness of the network strands 
m ~ Vh£, [S], as reported in Fig. [ija. 

With these experimental parameters we now can test 
the validity of our assumptions. Goo = Gs and Go = Gb- 
With g — 22.0 ± 3.5 we find that eq. 1 and 2 describe 
the volume fraction dependence of both moduli remark- 
ably well, as shown in Fig 3 where the lines correspond to 
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FIG. 3. Volume fraction dependence of high Goo (B) and low 

frequency modulus Go (O)- Lines are fits to the data using 
eq|T]and[2]with ga=77±5. 
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FIG. 4. Volume fraction dependence of the volume fraction 
occupied by the glassy network h and the strand dimensions 
(f, m). a) h increases linearly with (j), h = atp + b with 
a=3.5±0.3 and &=-0.15±0.03. b) f (■) is independent of (j), 
while m (Q) increases with </!> as m ~ (^°, with a=0.55±0.05. 



the theoretical prediction according eq 1 and 2. The re- 
markable agreement between prediction and experimen- 
tal data supports the idea that Goo is actually dominated 
by the contribution of the stretching modes, while Go is 
dominated by the bending modes of the network strands. 
However, we note that the value we obtain for the param- 
eter 5 ~ 22 is larger than what is typically estimated in 
theory [211 l2I]- This large value for g reflects the re- 
markable strong dependence of the moduli on (j). Such 
increased dependence in comparison to that theoretically 
expected for porous media is most likely due to the struc- 
tural differences between our very course networks and 
a porous medium comprising dilute and non-interacting 
pores, as assumed in our modeling. Moreover, we also 
neglected the heterogeneities inside the stands. Assum- 
ing that their contribution to the elastic modulus would 
linearly depend on </), this would lead to an enhanced de- 
pendence of Gs on (p and to a value g that is not only 
due to the pore, i.e. strands geometry. 

Nonetheless, the good agreement between the mechan- 
ical properties of the arrested spinodal network and the 
scaling behavior predicted by our simple model supports 
the condition that stretching and bending modes are de- 
coupled because of a delayed response of the bending 
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mode. This feature is difBcult to understand, but as 
denoted before could be related to the large differences 
in the length scales involved in the mescoscale organi- 
zation of the network. Indeed, the bending of the thick 
strands may entail structural rearrangements of a large 
number of particles, which is a slow process. In con- 
trast, stretching only involves local stretching of particle- 
particle bonds, which occurs instantaneously upon appli- 
cation of a stress. 

Note, that the difference between particle and strand 
size also dictates the unusually strong dependence of the 
elastic modulus on volume fraction. Because of this dif- 
ference the gels exhibit features of porous media, the 
elasticity increasing exponentially with instead of the 
typical power law increases observed in classical col- 
loidal gels. Indeed, for colloidal systems, where spin- 
odal decomposition is for instance induced by depletion 
H [3 ini [ini m HB HO] , the arrest of the phase separation 
process generally leads to gels where the ratio between 
particle and strand size barely exceeds one order of mag- 
nitude. Both, bending and stretching modes should here 
still contribute to the elasticity of the system; however, it 
is unlikely that both modes will be as clearly separated 
in time, as this is the case in the lysozyme gels. For 
other protein-gels that are formed by spinodal decom- 
position, it appears likely that the mechanical response 
function will exhibit features similar to those observed 
here. Indeed, studies exploring the phase behavior of a 
series of different proteins [TT| indicate that a large dif- 
ference between protein and strand size may be a gen- 
eral feature of protein-systems that are deeply quenched 
into the liquid-liquid coexistence range. If this is indeed 
the case the resulting mechanical properties should be 
similar to those of the deeply quenched lysozyme sys- 
tems. To make a definite assessment, however, further 
investigations of the volume fraction dependence of the 
structural and mechanical properties of deeply quenched 
protein systems are required. 

In conclusion, we have investigated the rheology of the 
complex network structure that forms as a result of an 
arrested spinodal decomposition in solutions of the globu- 
lar protein lysozyme. For these systems the characteristic 
correlation length of the network is found to be remark- 
ably large as compared to the particle size. The mechan- 
ical response are characterized by two distinct moduli. 
Go and Goo that are well separated in time. Go and Goo 
being respectively associated to the long and short time 
response of the system. Using a porous medium approx- 
imation, we develop a simple model that links the struc- 
tural features of the network to its rheological proper- 
ties and their 0-dependence. The remarkable agreement 
between theoretical prediction and experiment indicates 
that Go and Goo are dominated by respectively the bend- 
ing and the stretching modes of the arrested network. 
The two moduli are distinct and measurable because of 
the hierarchical mesoscale organization separating the 
network characteristic length {fim) and the protein size 
(nm): this length-scale difference appears responsible for 



the separation of the time scales over which bending and 
stretching occurs. Our combination of experiments and 
modeling creates a framework for understanding the re- 
lation between rheology and structure. Our results may 
provide new guidelines to design amorphous solid-like 
network with specifically defined and controlled charac- 
teristic structural and mechanical properties with poten- 
tial applications in materials and food science 22, 31 , 32j. 
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EPAPS: determination of Go and Goo using the 
scaling behavior of 
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FIG. 5. Scaling of the inverse compliance on a master 
curve J*~^. a) Examples of the time evolution of the inverse 
compliance J{t)~^ presented in Fig. 2 in the article. From 
top to bottom: T=10 °C, (/> = 0.171, 0.149, 0.13, 0.11, 0.093 
and 0.075. J^^{(f) ~ 0.11) is highlighted as it serve as the 
reference for the scaling procedure, b) Master curve, J*~^ 
based on the reference J~^{4> — 0.11), showing the scaled 
inverse compliance for the examples shown in a). 



As shown in the EPAPS Fig. la (Fig. 2 in the arti- 
cle), the inverse compliance data not always extend to 
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sufficiently long times to be able to properly determine 
the terminal modulus Go , which should be determined in 
the time range where J^^ become time independent. To 
also determine Go for the data-sets with a somewhat lim- 
ited time window and to improve the relative accuracy of 
Go and Goo obtained for the different samples, wc use a 
scaling approach. In this scaling approach the functional 
form of the material response function of all our samples 
is presumed to be the same, differing only in the parame- 
ter Go, Goo and Tc, where Tc characterizes the time scale 
describing the transitional range between Goo at short 
time scales and Go at long time scales. In praxis we fix 
the data set obtained with J~\{(j>) = J~^{4' = 0.11) as 



reference and then scale all other data sets onto this curve 
by applying the following algorithm to J~^: 

r-\<l>) = J-}j{c^) 

_ J-^W-Go(^) 

~ Goo(</.)-Go(<^) ^^"■'■^^ ~ ^"•'"^^J + 

For the time axis we use t = Tc/Tref- As shown in the 
EPAPS Fig. lb, the scaling procedure works remarkably 
well supporting the hypothesis that the functional form 
of the material response function is the same for all our 
samples. We therefore use this scaling approach to de- 
termine Go and Goo and the results are reported in Fig. 3 
in the article. 
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